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ABSTRACT 

The distribution of secondary star masses in present-day post-common en- 
velope binaries (PCEBs) is calculated using four different models for angular 
momentum loss (AML) during the post-CE phase: only gravitational radiation 
(GR), GR + disrupted magnetic braking (DMB), GR + reduced MB, and GR 
+ intermediate MB. For the DMB model, we find that the number of PCEBs 
decreases abruptly by 38% once MB begins to operate for non-fuUy convective 
secondaries. We do not find a similar feature in the distributions calculated us- 
ing any of the other three AML models in which MB is not disrupted. This 
percentage decrease in the number of present-day PCEBs predicted using the 
DMB model is easily large enough so that an observed distribution of secondary 
masses or even spectral types in PCEBs can provide an important test of whether 
magnetic braking is indeed substantially reduced in secondary stars that are fully 
convective. We discuss briefly the feasibility of such observations. 

Subject headings: binaries: close — stars: low-mass — stars: evolution — stars: mag- 
netic fields — white dwarfs 



1. Introduction 



The angular momentum loss (AML) mechanisms that operate to drive mass transfer in 
cataclysmic variables (CVs) have been the subject of debate recently. Historically, it has been 
believed that gravitational radiation (GR) is the main driver for CVs having orbital periods 
of 3 hours or less and magnetic braking (MB) the main driver for orbital periods of 3 hours 
or greater. In the orbital period distribution of observed CVs, there is a significant shortage 
of systems with periods between 2 and 3 hours. This feature is referred to as the "period 
gap." A generally accepted model for the period gap, proposed independently by Rappaport, 
Verbunt and Joss (1983; hereafter referred to as RVJ) and Spruit & Ritter (1983), is that 
MB is severely reduced once the secondary star becomes fully convective. This occurs for a 
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mass of ~ 0.3-0.4 Mq, which corresponds to an orbital period in CVs of ~ 3 hours. Once 

MB ceases, the secondary, whose radius is shghtly larger than its main sequence radius due 
to the high mass transfer rates caused by MB, relaxes back to its thermal equilibrium radius 
and becomes detached. Mass transfer recommences once GR brings the Roche lobe back into 
contact with the secondary at ~ 2 hours. This model has become known as the disrupted 
magnetic braking (DMB) model. 

Recently, the DMB model has been called into question by observations of magnetic 
activity in single stars in young open clusters (Sills et al. 2000; Pinsonneault, Andronov & 
Sills 2002). Sills ct al. (2000) find no significant decrease in magnetic activity through the 
mass range where the secondary becomes fully convective. They also find that the magnetic 
activity saturates at stellar rotation rates larger than some critical rotation rate, flcrit, a-nd 
claim that the prescription used in RVJ overestimates the strength of the AML due to MB at 
high rotation rates. Andronov et al. (2003) have constructed a MB prescription that better 
fits the observational data from open clusters. In their model, Jmb varies as (as in RVJ) 
for rotation rates slower than flcru, and as Q for faster rotation rates. This model, which 
has become known as "reduced magnetic braking (RMB)," if applicable to the secondary 
stars in CVs, poses serious problems for the DMB model. However, we should also note 
that the RMB model provides no alternate explanation of the period gap or of the observed 
brightnesses of CVs above the period gap. 

Ivanova & Taam (2003) have recently proposed a MB prescription based on a two- 
component coronal model of stellar magnetic fields (Mcstcl & Spruit 1987) and on observa- 
tions of X-ray emission from rotating dwarfs that provide an empirical relationship between 
stellar activity and rotation (Pizzolato et al. 2003). Their model is similar to the RMB 
model in that there is a different dependence upon Q for fast rotators than for slow rotators. 
However, the dependence of J on Q in their model for fast rotators is a bit stronger than 
in RMB (see section 2.2.). Hence, we shall designate this model as "intermediate magnetic 
braking (1MB)." 

The same AML mechanisms that drive mass transfer in CVs also act prior to the CV 
phase, when the binary is detached. Close binary stars in this stage of their evolution are 
typically referred to as "post-common envelope binaries (PCEBs)." To avoid any potential 
ambiguity, we note that in this paper the term PCEBs specifically refers to binaries that (1) 
are detached, (2) contain a white dwarf (WD) primary and a main-sequence or brown dwarf 
secondary, and (3) have undergone a single CE phase. We present in this paper population 
synthesis calculations of the secondary mass distribution in present-day PCEBs for the DMB, 
RMB and 1MB models, as well as for no MB (i.e., only GR). We show that this distribution 
can provide a relatively straightforward observational test of whether magnetic braking is 
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indeed disrupted once the secondary star becomes fully convective. 

We note that most early population synthesis calculations of PCEBs focused on explor- 
ing statistically the possibility that newly-formed PCEBs could be responsible for bipolar 
planetary nebulae (deKool 1990; deKool & Ritter 1993; Yungelson et al. 1993; Han et al. 
1995). With one exception, these studies did not model the subsequent evolution of the 
newly-formed PCEBs and thus did not consider AML mechanisms. Only deKool & Ritter 
(1993) calculated models of the present-day PCEB population using DMB and we compare 
our results with theirs in section 4. Recent investigations of PCEBs have been done by 
Schreiber & Gansicke (2003), Han et al. (2003), Willems & Kolb (2004), and Nelemans & 
Tout (2005), but none are suitable for comparison. Schreiber & Gansicke (2003) did not 
perform full population synthesis calculations, Willems & Kolb (2004) computed formation 
models of PCEBs, but not present-day, evolved models including an AML phase, Han et al. 
(2003) only considered sdB binaries, and Nelemans & Tout (2005) were mainly concerned 
with reconstructing the evolution of observed PCEBs through prior phases of mass transfer 
and assumed negligible orbital evolution of the PCEBs due to AML. 



2. Method 

2.1. Population Synthesis Code 

The Monte Carlo population synthesis code used in this study is the same as the one used 
in Politano (2004) to model ZACVs. This code is described in Politano (1988, 1996, 2004), 
and we must refer the reader to those papers for detailed discussions. Here we summarize 
only the key assumptions and features of the code and discuss how the code was used to 
calculate the present-day PCEB population. 

For all calculations, we begin with 10^ zero-age main sequence (ZAMS) binaries. Fol- 
lowing a standard approach, we assume that the distribution of these ZAMS binaries can 
be written as a product of three separate distributions over primary mass, mass ratio, and 
orbital period. We use a Miller & Scalo (1979) distribution for the primary masses, a distri- 
bution that is fiat in q for the mass ratios (i.e., g(q) dq = 1 dq, where q = M^/Mp, Duquennoy 
& Mayor 1991; Mazeh et al. 1992; Goldberg et al. 2003), and a distribution that is flat in 
log P for the orbital periods (Abt 1983). For a given primary mass, the secondary mass 
is chosen from the distribution, F(Ms) = f(Mp)g(q), where f(Mp) is the above-mentioned 
Miller & Scalo (1979) distribution used for the primary masses. 

To model the population of PCEBs, a given ZAMS binary is evolved to the point where 
the primary contacts its Roche lobe during its ascension of one of the giant branches. Wind 
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loss during ascent of the giant branch(es) is incorporated via a Reimers (1975) prescription. 
Relationships between the radius of the giant, its core mass and its total mass are derived 
from analytic fits to detailed stellar evolution models (Politano 1988, 1996). Simple energetic 
considerations are used to relate the pre- and post-CE orbital separations for the CE phase 
and a standard constant acE prescription is used with acE = ^ (see Tutukov & Yungelson 
1979; Politano 2004). The classical prescription for GR (Landau & Lifshitz 1951) and various 
prescriptions for MB (see next section) are used to describe AML during the post-CE phase. 
We use detailed stellar models from the Lyon group (Chabrier & Baraffe 1997; Baraffe et 
al. 1998, 2003; Chabrier et al. 2000) for low-mass secondaries (< 0.5 M©) and fits to stellar 
models from Webbink (see Politano 1988, 1996) for secondaries with masses greater than this. 
In our models, the secondary is fully convective for <0.37 Mq. To avoid any possible 
confusion with CV evolution, we emphasize that the secondary stars in PCEBs do not lose 
mass as a result of mass transfer since the binary is detached. Therefore, the secondaries in 
PCEBs do not "become" fully convective. They either are fully convective or are not fully 
convective in our models, depending on whether their mass is less than 0.37 Mq or greater 
than 0.37 Mq, respectively. Finally, we terminate MB in all AML models once the secondary 
develops a radiative envelope at 1.25 Mq. 

In calculating the population of present-day PCEBs, four timescales are of interest: (1) 
tb, the time that the progenitor binary was formed (measured from the beginning of the 
Galaxy), (2) tev,p, the time it takes the primary to evolve off of the main sequence, become 
a giant, and contact its Roche lobe to initiate the CE phase, (3) tpcEB, the time from the 
end of the CE phase until the present epoch, and (4), toai, the age of the Galaxy. These 
four timescales must satisfy the following constraint, th + tev,p + tpcEB = tcai- We make the 
following assumptions regarding these time scales: (1) the stellar birth rate throughout the 
Galaxy's history has been constant, (2) the CE phase is so rapid that the time spent in it 
by the binary is negligible compared to the other time scales (Meyer & Meyer-Hofmeister 
1979), and (3) the age of the Galaxy is 10^° yrs. 



We use four different AML prescriptions that differ in how MB is modeled. Model 1 is 
no MB, only GR. We use the standard prescription for GR (Landau & Lifshitz 1951): 



where G is the gravitational constant, c is the speed of hght, a is the orbital separation. Mi 
and M2 are the masses of the two stars, and Q is the orbital angular frequency. Model 2 is 



2.2. AML Prescriptions 




(1) 
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GR + RMB, where the prescription for RMB is given by Andronov et al. (2003): 



and where M and R are the mass and radius, respectively, of the secondary star and Qcru is 
the rotation rate at which MB saturates. Model 3 is GR + 1MB, where the prescription for 
1MB is given by (Ivanova & Taam 2003): 

0.5 
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In this prescription, fl^ is the rotation rate at which the dependence of the X-ray luminosity 
on D, changes form, which we assume is equal to 10 Qq, and we assume is equal to Tfj © 
(Ivanova & Taam 2003; R. Taam, priv. comm.). Finally, model 4 is GR + DMB, where the 
prescription for DMB given by RVJ is used: 

"^^^^"lo, if M< 0.37 Mo. 



3. Results 

In Figure 1, we show the distribution of secondary masses in present-day PCEBs for 
each of the four AML models described in section 2.2. The dotted hue is the distribution 
for GR only, the dashed hue for GR + RMB, the dashed-dotted line for GR + 1MB, and 
the solid line for GR + DMB. The scale on the y-axis has been arbitrarily normalized to 
facilitate the discussion below. 

The most striking feature in this figure is the sharp drop in the number of present-day 
PCEBs once MB is turned on at M^ = 0.37 M© in the DMB model. The distributions 
calculated using the other AML prescriptions do not show a similar feature. This sharp drop 
is due to the increased efficiency of MB to bring systems into contact compared to GR in 
the DMB model. Rewriting equation 1 in units similar to equation 4, we have 

4„ = -a.2xI0^».,„.™(j,J^)^(i^)'^ (5) 

Comparison of equations 4 and 5 shows that J is ~ 50 times greater for MB than for GR 
in the DMB model (assuming a typical WD mass, Mi = 0.6 Mq). As a result, when MB 
is operating, PCEBs are brought into contact relatively quickly following the CE phase and 
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become CVs before the present epoch is reached. On the other hand, if only GR is operating, 
many more PCEBs remain detached at the present epoch. 

The percentage decrease in the number of PCEBs once MB is turned on in the DMB 
model is 38%. This percentage decrease compares quite favorably with that found by deKool 
& Ritter (1993) (~44%, see their Fig. 4c), who also calculated the secondary mass distri- 
bution in present-day PCEBs using DMB. 

For the DMB model, we also investigated the effects of different values of acE and 
different ZAMS mass ratio distributions on the resulting present-day PCEB secondary mass 
distribution. In addition to acE — 1 and g(q) = 1, we calculated distributions for acE = 
0.6, 0.3 and 0.05, and for g(q) oc q and g(q) cx q~°'^. We find that while the total number of 
PCEBs is sensitive to the choice of acE^ the percentage decrease in the number of PCEBs 
at 0.37 Mq remains relatively unaffected for acE between 1.0 and 0.6, and then increases 
from 37% to 73% as acE is reduced to 0.05. We also find that while the overall shape of 
the distribution is affected by the choice of g(q), the sharp drop at M^ = 0.37 Mq remains 
prominent in each distribution, with only a slight variation in the percentage decrease of 
PCEBs (29% to 37%). 

4. Discussion 

In our population synthesis calculations using the DMB model, MB does not operate 
for secondary stars that are fully convective (M^ < 0.37 M© in our models). For secondary 
stars that are not fully convective (0.37 M© < M^ < 1.25 Mq in our models), MB does op- 
erate in the DMB model. The number of present-day PCEBs predicted using the DMB 
model decreases abruptly by 38% when magnetic braking is turned on at M^ = 0.37 Mq. 
The present-day PCEB secondary mass distributions predicted using the other three AML 
models, in which MB is either always off (GR only) or always on (RMB and 1MB) show no 
similar abrupt decrease. Notably, the percentage decrease predicted using the DMB model 
is insensitive to the assumed choices for the common envelope efficiency parameter and the 
distribution of mass ratios in ZAMS binaries, two input factors whose uncertainties have 
notoriously plagued population synthesis calculations of close binary systems. Such insensi- 
tivity implies that this prediction is fairly robust. Consequently, observations of PCEBs with 
secondary masses in the range ~ 0.25 - 0.5 Mq (approximately M and late K spectral types) 
should provide an important test of whether MB is indeed disrupted once the secondary 
becomes fully convective. 

As few as six years ago, there were only ~ 40 known PCEBs (Hillwig et al. 2000) and 
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the prospects for observational comparison with theoretical models might have seemed bleak. 
However, with the advent of large observational surveys such as the Sloan Digital Sky Sur- 
vey (SDSS), the observed number of detached WD/MS binaries has increased substantially. 
Silvestri et al. (2006) have recently extended the work of Raymond et al. (2003) and have 
compiled a catalog of 746 spectroscopically-identified, detached close binary systems contain- 
ing a WD and a low-mass MS secondary star in the SDSS through the Fourth Data Release 
(Adelman- McCarthy et al. 2006). They find that the vast majority of these systems (~ 700) 
contain M dwarf secondaries and that the distribution of secondary spectral types in these 
binaries peaks at a spectral type of M4 and trails off quickly for earlier and later spectral 
types (see Fig. 10a in Silvestri et al. 2006). However, the majority (~ 80%) of the detached 
WD/MS binaries in the Silvestri et al. (2006) sample were targeted for spectroscopy because 
their colors resembled other, higher priority objects, such as quasars. Richards et al. (2002) 
discuss in detail how quasar candidates were selected in the SDSS and the regions in color 
space that were excluded because of potential contamination of the quasar sample. One such 
excluded region is precisely where WD-M dwarf binaries reside. Consequently, as noted in 
Silvestri et al. (2006), their sample is neither well-defined photometrically nor statistically 
complete and should not be considered as representative of the secondary mass function in 
PCEBs. 

Nevertheless, this recent large increase in the number of detached WD/MS binaries is 
encouraging and suggests that the test of the DMB model we propose is feasible obervation- 
ally. This feasibility is greatly aided by the following factors: (1) the predicted percentage 
decrease in the number of PCEBs in the DMB model is large enough so that a modest sample 
(~ 50-100 PCEBs) should suffice; (2) the range of secondary masses to be targeted is narrow 
and fairly well-defined (~ 0.25-0.50 M©); (3) precise determinations of the secondary masses 
are not necessary, spectral types will suffice, and (4) precise orbital periods arc not needed, 
only a coarse determination of whether the system could have undergone a CE phase (i.e., 
has an orbital period less than ~ 10 days). 

An observational challenge that will need to be addressed is the detection of PCEBs 
containing either a hot WD and a very low mass secondary or a cool WD and an early M /late 
K secondary. In the former case, the WD is much brighter than the secondary in the optical 
and the system may be mistaken for a single WD based on its colors. In the latter case, the 
WD is hidden by the glare of the secondary, and the system may be misidentifed based on its 
colors as a single early M or K star. For example, Schreiber & Gansicke (2003) computed the 
colors expected for PCEBs containing a WD in the temperature range 6000 - 35,000 K and a 
ZAMS secondary with spectral type KO to M6. They compared their simulated colors with 
the U - B selection criteria for the Palomar-Green survey (Green et al. 1986) and concluded 
that PCEBs containing a relatively cool (Tg// < 15,000 K) WD would have been included 
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in the Palomar- Green survey only if they had a companion of spectral type ~ M4 or later. 
We note that observations in the UV or blue region of the optical spectrum (~ 3500-5000 
A), where early M dwarfs contribute very little flux, may help to reveal a WD companion if 
it exists. However, optical observations in the blue would need to be at a higher resolution 
than SDSS to achieve an acceptable signal to noise ratio (N. Silvestri, priv. comm.). EUVE 
or X-ray observations may prove useful in revealing K stars that are hiding a WD companion 
(e.g., Green, Ah, & Napiwotski 2000; Good et al. 2005). 

Warm thanks to Drs. C. Garden, D. Drapes, D. Hoard, N. Silvestri, K. Simkunas, P. 
Szkody, R. Taam and B. Willems for very useful discussions about this work. We are partic- 
ularly grateful to Dr. N. Silvestri for providing access to her SDSS observations in advance 
of publication. Finally, we thank the referee for very helpful comments that improved the 
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Fig. 1. — Theoretical distributions of the secondary mass in present-day PCEBs for the four 
different assumed AML models described in the text: GR only (dotted line), GR + RMB 
(dashed Une), GR + 1MB (dashed-dotted hne), and GR + DMB (solid line). The y-axis has 
been arbitrarily normalized to facilitate comparison. Corresponding spectral types are shown 
along the top axis. Spectral type-mass relationships from Kirkpatrick, Henry, & McCarthy 
(1991) and Kirkpatrick & McCarthy (1994) were used for M dwarfs. 



